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Abstract In the present article we intend to accomplish a rigorous analysis which allows to 
elaborate useful symbolic and numerical codes leading to an accurate evaluation of the thermal 
radiation. We will consider the case of the electromagnetic field in axially symmetric systems, a 
case of real interest in cryogenics. 



PACS numbers: 41.20.-q 

I. ELECTROMAGNETIC FIELD IN CONDUCTIVE BOUNDED CYLINDRICAL SYSTEMS 

The first topic is the electromagnetic field inside a hollow cylinder cavity with conducting walls and the spectrum 
of resonance frequencies associated to the various modes of the field. 

We will assume that the field has a monochromatic time dependence e*"^*. The symmetry axis will be taken as the 
z-axis of the coordinate system. With this choice, the field will split in two orthogonal components as 

E = e,E, + E± (1.1) 
The walls of the cavity are situated at p = b, z — and z — I. 

The tangential component of the field should vanish on the conducting walls, that is Ez = E^ = on the cylinder 
surface p = h and Ep = E^ = Q aX z = Oandz = I (fig.l). Inside the hollow cavity the field will obey the wave equation 
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or 
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S7' + -Je = Q 

for the case of harmonic field of frequency uj. The axial component of the space part of the electric field obeys the 
same equation: 



V' + —jE,ip,cP,z)^0 (1.2) 
and so does the transverse part E± : [Ep, Eff,). The equation can be written in cylindrical coordinates as 

f^ + lA + l^ + ^ + ^^ij ^0 (13) 
ydp^ pdp 90^ dz'^ J 

and assumes a separate solution of the type 

E,{p,^,z)=R{p)<^{^)Z{z) 

. Thus, the Helmholtz PDE above becomes 
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R \ dp"^ P dp J p"^^ dcj)^ Z dz'^ 



(1.4) 



The reflections at the top z = /and bottom z = surfaces impose a z-dependence appropriate for standing waves. 
Thus, in the eq. (1.6) the constant fc^ > and the solution Z{z) will be a combination of sinfczandcosfcz which 
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will obey the boundary conditions. Further, in order to have a uniform (j) dependence, $(0) should be of the form 
$(,/,) = e^*""^ with integer m = 0, 1, 2, 3, ... 
Thus, the radial equation becomes 

(fR IdR f ^ _^ ^ 

+ [^^ - — \r = Q (1.5) 



dp^ p dp \ p^ 

where 7^ = ^—k'^. The equation (1.7) is the Bessel equation of order m = integer and the solution is a combination 
of Bessel and Neumann functions of order m , Jm{'yp)a,ndNm{'Jp)- If we take into account that the solution should 
be finite for any p, including p = and remember the behavior of Bessel and Neumann functions in the vicinity of 
the origin 



•^'"^''^^ r(m + i) (2) 



r(m + i) 

respectively 
for m > 1 and 
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7Vo(a;)^- [ln(|)+0.5772...' 



than only Jmiip) is acceptable, because the second solution is singular at the origin. In the case of Transverse 
Magnetic (TM) modes =0 everywhere in the cavity, while the transverse components of the electric and magnetic 
field are 



and, respectively 



E. = -V^^ (1.6) 
7 oz 



B± = X V^i?. (1.7) 



where V_l : is the Hamilton's operator in polar coordinates. 

Than must assume the form 

E,{p, <j>, z) = ^0 J„(7p)e±^'"^ cos (p7r|) withp = 0, 1, 2, ... (1.8) 

in order that E± to satisfy the boundary conditions at z = Oandz = I. Also, the boundary condition = for 
p = b requires that the argument 76 should be one of the zeros of the Bessel function J^- Thus, 7 will assume only a 
set of discrete values 
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where Xmn is the n-th zero of Jm- 

The general solution of the field equations (1.10) inside the cavity that obeys the appropriate boundary conditions 
will be a superposition of the normal modes : 
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E,{p,^,z) = ^ ^^i;„npJm(a:™„^)e±''"^cos (pirj) (1.9) 

m=0 n=l p=0 

The resonance frequencies of the cavity will depend on three parameters and will be given in the case of TM modes 

by 



with m=0,l,2,... n=l,2,... p=0,l,2,... 



FIG. 1: The geometric parameters of the system 



II. HOLLOW CYLINDRIC RING WITH CONDUCTING WALLS 



The second topic is the electromagnetic field inside a hollow cavity of the form of a cylindric ring with conducting 
walls, with inner radius a, outer radius 6and height I. 

The PDE and the boundary conditions that the fields obey will be the same as in the case of the cylinder cavity, 
except that a supplementary boundary at p = a{> 0) appears. Since the origin p = no longer belongs to the domain 
of values that p assumes, than the particular solution of the radial equation will be now a linear combination of both 
Bessel and Neumann functions: 



In order that the field to satisfy homogenous boundary conditions on both p — aand p — b boundaries, the radial 
funtion above should assume the form 



Rrnip) = AmJmi'Jp) + BmNmilp) 



(2.1) 



R{p) = A Jra{lmnp) 




NmilninP) 



(2.2) 



where "fmnji^ = Ij 2, 3, ... are the solutions of the equation 



They will be determined by appropriate numerical methods. 

In the case of TM modes, the general solution for the axial component is the linear combination: 
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E^{p, ^) = V V V Emnp JmilmnP) " , \ Nm{lmnP) 6^"""* COS (pVT- 



(2.3) 



while the resonant frequencies are 




(2.4) 



with m=0,l,2,... n=l,2,... p=0,l,2,. 
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III. CONCLUSIONS 

The presented formulae allow an exact calculus of the electric field inside axial cavities with conducting walls, 
which may be used for a rigorous evaluation of the radiation thermal flux in such a system. For cryogenic systems, the 
numerical calculations need such accurate formulae, taking into account that a large number of coaxial insulators will 
be used and any errors will propagate and will be recursively amplified. The only way to minimize the global errors is 
to keep them as low as possible at any stage, thus involving from the beginning these exact analytical solutions. That 
is why the numerical calculus of the special functions involved must be carefully conducted, implying an increased 
computing effort in order to precisely describe the total radiation heat transfer rate in cryogenic installations. 



APPENDIX 



The Bessel functions J^{x)aYe solutions of Bessel equation 
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dx'^ X dx 
and assume the following series expansions 
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R{x) = 
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(A.2) 
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J=0 



{-ly (x\^i 



1) 



They satisfy the following orthogonality relation on the interval p G [0, a] 
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(A.3) 



The above solutions are independent for all u except the case of integer = m = 0, 1, 2... when 

J_„(X) = i-irJmix) 

In this case a second independent solution of Bessel equation is the Neumann function 



N^x) = 



Jv{x) COS V1T — J-v{x) 



sm i/TT 

Both Bessel and Neumann functions as well as their combinations 
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Hl^\x) = J,{x)+iN,{x) 



and 



H^^\x) = J,{x)-iN,{x) 
i.e the Hankel functions, satisfy the recursion formulae 



2v 

X^^i{x) + X^+i{x) = —X^{x) 

X 



and 



(A.5) 



X,_i(,xO =2 
where X(x) denotes a Bessel, Neumann or Hankel function. 



dX,{x) 
dx 



(A.6) 
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